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Abstract 

Interquark confinement potential is calculated in the dual Monopole Nambu- 
Jona-Lasinio model with dual Dirac strings suggested in Refs.[l,2] as a functional 
of a dual Dirac string length. The calculation is carried out by the explicit inte- 
gration over quantum fluctuations of a dual-vector field (monopole-antimonopole 
collective excitation) around the Abrikosov flux line and string shape fluctuations. 
The contribution of the scalar field (monopole-antimonopole collective excitation) 
exchange is taken into account in the tree approximation due to the London limit 
regime. 
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1 Introduction 



The dual Monopole Nambu-Jona-Lasinio model (MNJL) with dual Dirac strings as con- 
tinuum space-time analogy of Compact Quantum Electrodynamics (CQED) [1] has been 
formulated in Ref. [2-4]. The MNJL-model is based on a Lagrangian, invariant under 
magnetic U(l) symmetry, with massless magnetic monopoles self-coupled through a local 
four-monopole interaction [2,3]: 

£(x) = x(?)i'fdvX(B) + G[x{x)x{x)] 2 - Gilxix^^xi^MXixh^xix)}, (1-1) 

where x{ x ) is a massless magnetic monopole field, G and Gi are positive phenomenolog- 
ical constants. Below we show that G± = G/A. The magnetic monopole condensation 
accompanies itself the creation of massive magnetic monopoles Xm{x) with mass M, Xx~ 
collective excitations with quantum numbers of the scalar Higgs field a with the mass 
M a = 2M and the massive dual-vector field with the mass Mq defined as [2,3]: 

2 2 

M2 ° = 2G[ ~ h^ [Jl{M) + M2j ^ M ^ ^ 
where J\{M) and ^(M) are quadratically and logarithmically divergent integrals [1,2] 

J 'M = /Sm^- a2 -^"( 1 + ^)' 

Here A is the ultra-violet cut-off. The mass of the massive magnetic monopole field 
Xm{x) obeys the gap-equation [2,3]: 



M = -2G< x(0)x(0) >= ^Ji(M) (1.4) 



derived from the effective Lagrangian of the scalar cr(x) and the dual-vector C^{x) fields 
by virtue of the suppression of the direct transitions a < — > vacuum. On the other 
hand, due to one-loop corrections to the mass of the monopole field derived by using the 
Lagrangian Eq.(|l.l|) the gap-equation should read 

M = -2^G + G^j< x(0) X (0) > . (1.5) 

Since the level of the collective xx-excitations should be completely compatible with the 
monopole level, the gap-equations Eq.( [l~2|) and Eq. (|1.4|) should coincide. This fixes G\ 
in terms of G as G\ = G/A. 

As has been shown in Refs.[2, 3] the vacuum expectation values of time ordered prod- 
ucts of densities expressed in terms of the massless monopole field, i.e., the magnetic 
monopole Green functions 

G(x u ...,x n ) =< 0|T(x(a;i)r 1 x(xi) . . -xM F n X > CO nn. , (1-6) 
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where r,(z = 1, . . . , n) are the Dirac matrices, are given by [2,3] 

G(x 1 , ...,x n )=< 0|T(x(xi)rix(xi) . . .x(x n )T n x(x n ))\0 > conn .= 
= (M) < 0|T(xAf(a;i)riXM(^i) • • • XmMTuXmM 

x expij d 4 x{-gxM{x)l v XM{x)C v {x) - kxm{x)xm{x)o{x) 

+£ int [a(x)]})|0>W. • (1-7) 



Here |0 > (M) is the wave function of the non-perturbative vacuum of the MNJL-model 
in the condensed phase and |0 > the wave function of the perturbative vacuum of the 
non-condensed phase. Then, £ int [cr(a;)] describes self-interactions of the cr-field: 

AntKx)] = -KM a a\x) - ^cr 4 (x). (1.8) 

The self-interactions £i nt [cr(x)] provide cr-field loop contributions and can be dropped out 
in the tree cr-field approximation accepted in Refs. [2-4] . The tree cr-field approximation 
can be justified keeping massive magnetic monopoles very heavy, i.e. M 3> Mc- This 
corresponds to the London limit M CT = 2M > Mc in the dual Higgs model with dual 
Dirac strings [5-7]. The inequality M a ^> Mc means also that in the MNJL model we 
deal with Dual Superconductivity of type I [4]. In the tree cr-field approximation the r.h.s. 
of Eq. ( |1.7| ) acquires the form 

G(X ± , ...,X n )=< 0\T(x(xi)Tlx(xi) . . . x(x n )T n x(x n ))\0 > C onn.= 

= (M) < 0|T(xA/(a;i)riXM(xi) . . . XM{x n )Y n XM{x n ) 

expi / d A x[ - gxM(x)YxM(x)C v (x) - kxm{x)xm{x)(t{x)})\$ >^ n . • (1-9) 



For the subsequent investigation it is convenient to represent the r.h.s. of Eq.( |1.9| ) in 
terms of the generating functional of the monopole Green functions [2-4] 



n 5 5 

G(Xl, . . . , X n ) = [] T~7 \^\f]i V] 



(1.10) 

?7=?7=0 



Y \f, v {x)F^{x) 



where 77(77) are the external sources of the massive monopole (antimonopole) fields, and 
Z[f],fj] is the generating functional of the monopole Green functions defined by 

Z [V,V) = 1^ J VxuVxhfDC^Vo- expz J d 4 x 

+ 1 - M 2 C C,(x) C»{x) + \ d,a(x) d^a(x) - I M 2 a a 2 (x) 
+XM(x)(iY <9 M - M - g^C^x) - kct(x)) Xm(x) 

+V(x) Xm{x) + Xm{x) T]{x) + £ frcc quark(^)J • (1.11) 

The normalization factor Z is defined by the condition Z[0, 0] = 1. The coupling constants 
g and k are related by the constraint 

2 2 

U J 2 {M) = ^-J 2 {M) = l (1.12) 



127T 2 ~ * v ~ ' 87T 2 
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or k? = 2g 2 /3 [2,3]. Then, £ free qU ark(^) is the kinetic term for the quark and antiquark 
U ee quark (x) = - £ m> / dT ( dX dr T) ^dT^ A ' 6(A){X " Xi(T)) - (L13) 



i=q,q 



In our consideration quarks and antiquarks are classical point-like particles with masses 
m q = m>q = m , electric charges Q q = —Qq = Q, and trajectories X^{r) and X^(t), 
respectively. The field strength F» v {x) is defined [2-4] as F» v {x) = £» u {x) - *dC» u (x), 
where dC^(x) = d»C v (x) - d v C»(x), and *dC u (x) is the dual version, i.e., *dC^(x) = 
^e^ va,3 dC a p(x) (e 0123 = 1). The dual "chromo" -electric field strength £^ u (x) ) induced by 
a dual Dirac string, is defined following [2-4] as 

where = X M (r, a) represents the position of a point on the world sheet swept by the 
string. The sheet is parameterized by internal coordinates — oo < r < oo and < a < n, 
so that X M (r, 0) = X^q{t) and X M (r, n) = Xq{t) represent the world lines of an anti- 
quark and a quark [2-7]. Within the definition Eq.( |1.14|) the tensor field £ tll/ {x) satisfies 
identically the equation of motion, d fl F lw {x') = J u (x). The electric quark current J u (x) 
is defined as 

r( x ) =Y.Q l j dr^j^-S^ix-X^r)). (1.15) 

i=q,q T 

Hence, the inclusion of a dual Dirac string in terms of £^ v [x) defined by Eq. (|1.14|) satisfies 
completely the dual electric Gauss law of Dirac's extension of Maxwell' s electrodynamics. 

The ground state of the massive dual-vector field C M (x) coupled to a dual Dirac string 
acquires the shape of the Abrikosov flux line [2-7] 

C[e(x)] = - 1 d 4 x' A{x - x 1 ) d*£» v (x' ), (1.16) 

where A(x — x') is the Green function 

/rl^h p * (>X X \ 

Wt Ml - J? - ,0 • W 

Integrating out the dual-vector field fluctuations c^(x) around the shape of the Abrikosov 
flux line, C M (x) = C tl [£(x)] + c^(x), and the scalar a-field [4] we obtain the generating 
functional of the monopole Green functions in the following form: 

z [v,v] = \\ v Xm^Xm expi j d 4 x \C^{xm{x), Xm(x), C v [£ (x)]} 

+Xm(x) {i^d^-M-gY C^[£{x)\) X m(x) + fj(x) X m(x) 

+Xm{x) ri(x) + £ free quark (a:) J , (1-18) 

where C cS {xm{x),Xm(x),C u [£(x)}} reads 

C cS {xm{.x),xm{.x),C v [£(x)}} = C st rm e {C' / [£(x)}} 

g 2 _ k 2 _ 

[Xm{x)i^Xm{x)} [xm{x)YXm{x)} + 77V7^[Xm(x)xm(x)] 2 . (1.19) 



2M 2LAMV 7,MAMV 7JLAMV 7 ' v /j ■ 2M 2 
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The Lagrangian of the dual Dirac string C stTing {C u [£(x)}} is defined [3-6] 

J d A xC stTing {C"[£{x)}} = \m 2 c J J d 4 xd 4 y£, a (x)A:(x-y,Mc)£^(y), (1.20) 

where A^{x - y, M c ) = « + 2d a d u /M 2 )A{x - y; M c ). 

The effective Lagrangian Eq. ( |1.19|) integrated over the massive monopole fields Xm(%) 
and Xm(x) defines the string energy, i.e. the interquark potential, as a functional of the 
string shape. 



2 Confinement potential 

The interquark confinement potential is related to the energy of the string which is defined 
as follows [2-7]: 

W = -J d 3 x£ stling {C»[£(x)]} + J d 3 x^ < 0|T(( - J^[x m (x)^ X m(x)} 



K 2 

x [Xm(x)^xm(x)] + ^[xm(s)xm(x)] : 



x exp-ig f d 4 y X M(y)rC,[£(y)}xM(y))\0 > (M) . (2.1) 



The interaction caused by the integration over the a-field fluctuations gives a trivial 
constant contribution to the energy of the string [4] and can be dropped out. In the 
momentum representation of the vacuum expectation values the energy of the string is 
then defined by [4]: 

w = -/ ** W*)l) - /A^ 

xJ _U 7 4/-% t X }. ( 2 . 2) 

1 M -h + gC[£(x)\ I J (2w)H \ " M -k 2 + gC[£(x)] J 

The momentum integrals have been calculated in Ref.[4]. This yields the energy of the 
string: 

W = -Jd 3 xjr, stIing {C[£(x)}} 

1 1 ( 9 " '■[J 1 (M)+M 2 J 2 (M)]\ J d 3 xC,[£{x)}C^[£{x)}. (2.3) 



2 Ml V8tt 21 



By using Eqs.([L2|) - ( |Q|) , the relations Gi = G/A and M a = 2 M we bring up the 
coefficient of the second term to the form 

2 2 — 

- ^[Ji(M) + M 2 J 2 (M)} =M 2 ~^- = M 2 + 8 g 2 ^—. (2.4) 

Thus, the energy of the string containing quantum fluctuations of the scalar and dual- 
vector fields around the shape of the Abrikosov flux line is given by 

W = -j d 3 xC string {C[£(x)]} 

~\ M2 °{ l + w ^r) 2 / d ' x ° Mx)] c " [£{x)] - (2 - 5) 



The computation of the r.h.s. of Eq. (|2.5| ) we perform for the static straight string of the 
length L directed along the z-axis. In this case the electric field strength £^ v {x) does not 
depend on time and is given by [6] 



S(x) = e z QS(x)5(y) 



0U--L 



z + -L 



(2.6) 



where a quark and an antiquark are placed at Xq = (0, 0, \ L) and X_q = (0, 0, — ~ L). 
The unit vector e z is directed along the z-axis and 8(z) is the step-function. The field 
strength Eq. (|2.6|) induces the dual-vector potential 



< C(x) >= -iQ 



d 3 k k x 



sin 



k z L \ ifc . x 



4vr 3 k z M^ + k 2 
For the static straight string the term — / d 3 x£ s t r mg{C 1/ [£(x)}} reads 

oo 

-J d 3 xC stfmg {C u [£(x)]} = --M 2 c J d 3 x J d 3 x' J dx' 



(2.7) 



£oi(z) \9i + 



d 2 



1 ' Mq dxidxi 



A(x - x' , x - x', M c ) S 0j (x' 



+E i0 (x)lg° + 



d 2 



Ml dx Q dx° 



A(x - x' , x - x\ M c ) S l0 (x' ) 



M c I d 3 x I d 3 x' Six) -S(x') ( 1 



M 2 dz\ 



A(x — x', Mi 



c) 



L/2 L/2 

l -Q 2 M 2 c j dz j dz> J 

-L/2 -L/2 



Q 2 M 2 



c 



AtT 3 

Q 2 M 2 



An 2 



dk z 
If 



dk z 
If 



sin 



sin 



k_X 

'kX 



dk z [l + 



1 + 



Ml, 



d 2 k, 



Jk z (z - z' ) 



(2tt) 3 Ml + k 2 + k 2 



1 + 



Ml 



Ml, 



d 2 k< 



Ml + k 2 + k 2 



dk 2 



Ml + ki + k 2 ' 



(2.8) 



where Aj_ is the cut-off in the plane perpendicular to the world-sheet swept by the string 
[2-7]. We identify A± with the mass of the scalar field, i.e., A ± = M a = 2M [2-7]. 
For a sufficiently long string we can integrate over k\ and get 



- Jd 3 xC stTing {C[S(x)]} 



Q 2 Ml 



ATI 2 



dk, 



sin 



x 



£n[l + 



MT 
Ml, 



£n\ 1 + 



1 + 



Ml, 



Ml n 



(2.9) 



where we have neglected k z relative to Aj_. Dropping the infinite constant contributions 
independent of L we obtain [6]: 

r O 2 M 2 T / M 2 \ 

- j d 3 x£ string {C»[£(x)]} = L^-^ inll + -^A +2E 1 (M C L) 
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M C L 



1-e 



-M C L 



An 



(2.10) 



where E\{MqV) is the Exponential Integral function. For the calculation of the integral 
over k z we have used the auxiliary integral 



dx 



sin 2 x 



x- 



£n[a 2 + 



x 



2n £na + 



7T 



1 -e 



— 2aa 



aa 



-2nE 1 (2aa), (2.11) 



where in Eq.( |2.1UD we have set a = 1 and a = McL/2. 

The first term proportional to L gives the string tension <To calculated in the tree- 
approximation [5]: 



0o 



Q 2 M ( 



c 



£n[l + 



in 



Ml, 



(2.12) 



The last term in Eq. ( |2.5| ) induced by the quantum fluctuations of the dual-vector field 
around the shape of the Abrikosov flux line can be reduced to the form [6]: 



--M 2 
2 



c 



8g 2 < xx > \ 2 

M 2 



Q 2 M ( 



c 



4tt 2 



-in 1 + 



1 + 



ig z <xx> 



M 2 M a 



d^xC^Six^C^Six)} 

2 oo 



dk z 
'k 2 



sm 



k z L 



ln\ 1 + 



Ml, 



Ml j, 



Q 2 M 2 



c 



*>g 2 < xx >V 



in 



tn\ 1 + 



MT 
M 2 C/ 



+ 2E 1 (M C L) - 



Ml 



M C L 



-M C L 



(2.13) 



Collecting the pieces together we obtain the energy of the dual Dirac string, the interquark 
potential, as a function of the the length of the string L: 



W 



Q 2 M ( 



c 



An 

+2E 1 (M C L) 



1 



W < XX > , 32c/ 4 < XX > ; 



M C L 



Ml 



1 



_ e -M G L 



Ml 



Ml J 
Q 2 e-McL 
An L 



ln\ 1 + 



MT 
Ml, 



(2.14) 



The term proportional to L describes a linearly rising interquark potential leading to the 
quark confinement and gives the expression for the string tension 



Q 2 M ( 



c 



An 



Ml 



M„ 



Ml 



Ml 



MT 
Ml, 



(2.15) 



The last term in Eq. (|2.14|) is the Yukawa potential. 

Matching the string tension Eq. (|2.15|) with the string tension o"o calculated in the tree- 
approximation Eq.( [2.12D we accentuate a tangible contribution of quantum fluctuations 
of the dual-vector field around the shape of the Abrikosov flux line. This agrees with 
the result obtained in Ref.[6] in the dual Higgs model. 
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3 String shape fluctuations 



The string shape fluctuations we define as usually [8,7] by X M — > X M + ^(X), where 
rj^(X) describes fluctuations around the fixed surface S swept by the shape T and obeys 
the constraint r] tJ, (X)\Qs = [8,7] at the boundary dS of the surface S. The integration 
over the //-field we perform around the shape of the static straight string with the length L 
tracing out the rectangular surface S with the time-side T [8,7]. Allowing only fluctuations 
in the plane perpendicular to the string world-sheet and setting r) t (t,z) = r] z (t,z) = 
[8,7], we arrive at the fluctuation action 5S^[r] x ,riy] [7,4] 

„ _ T/2 L/2 

3Q 2 A 2 r r 

5S N [r] x ,T] y ] = ^-i J dt J dz[r] x (t,z)(-A)r] x (t,z) + (x^y)}, (3.1) 

-T/2 -L/2 

coming from the term / d 4 x C s t T mg{C" [£ (x)}} defined by Eq.( |1.20|) , where A is the Laplace 
operator in 2-dimensional space-time 

d 2 d 2 

The term / d 4 x C fM [£(x)} C fJ, [£(x)] in Eq. (|2.5|) , induced by the quantum fluctuations of the 
dual-vector C M and scalar a fields around the shape of the Abrikosov flux line, does not 
contribute to the fluctuation action for the case of the static straight string. In order to 
show this we use the expression obtained in Ref.[4]: 

6CJ£(x)]C»[£(x)] = 



^2 



f f d 3 k d 3 q k x q x + k y q y . f k z L\ . ( q z L\ 1 



sin — — sin 



47r 3 47r 3 k z q z \ 2 J \ 2 J M 2 +k 2 

i(k + q) ■ x [{k x + q x )r] x (t, z) + (ky + q y )r] y (t,z)} _ ^ 



Ml + q 2 

The contribution to the fluctuation action is given by 
d 4 x5C^[£(x)]C^[£(x)} = 

n2 [.4 ffd 3 kd 3 qk x q x + k y q y , ( k z L\ . ( q z L\ 1 
Q I d x / / — — — — sm sin 1 1 



4vr 3 47r 3 k z q z \ 2 ) \ 2 J M 2 + k 2 

[(k x + q x )Vx{t, z) + (k v + q y )v y (t, z)\ _ A ( 3 4 ) 



1 i(k + q) ■ x ( i [(k x + q x )r) x (t, z) + (k y + q y )r) y (t, z)] 
M^ + q 2 



Integrating over x and y we get 
d i x5C^ 1 [£(x)}C^[£(x)} = 

Q 2 I dt [ dz [f*±*± + ky% sin ( ^) sin < q ' L 



-T/2 -L/2 



2n 2 2n 2 k z q z \ 2 J \ 2 J M 2 + k 2 



x -^-^ e* ^ + »)* " ^ + S (k x + q x ) 5(k y + q y ) 
Mq + q 

x (J [(h + Qx)Vx(t, z) + (k y + q y )v y (t, z)] _ A = Q ^ ^ 
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Thus, Eg. defines completely the fluctuation action induced by string shape fluctu- 
ations around a static straight string with length L. As has been shown in Ref. [7], the 
fluctuation action Eq.( |3.1|) gives a Coulomb-like universal contribution [8] to the energy 
of the string: 

W^string— shape = ^ , (3-6) 

where a s tri ng = vr/12 and a str ing = tt/3 for openned and closed strings, respectively. 



4 Conclusion 

We have shown that in the MNJL model with dual Dirac strings the quantum fluctuations 
of a dual-vector field C M and a scalar field a around the shape of the Abrikosov flux line 
give the interquark confinement potential in the following form 



f _ Q 2 M 2 ( 8g 2 < XX > , 32/ < X X> 2 ' 
tot 4tt [ Ml M a + Ml 



n I M 

£n\ 1 + 



Ml, 



+ 2E l{ M c L)- — U-e-^ 



Q & ^ ^string 

4^ L ~' 



(4.1) 



where a st rmg = n/12 and a str i ng = n/3 for openned and closed strings, respectively. This 
interquark potential resembles the result obtained in the dual Higgs model with dual 
Dirac strings [5-7]. Unlike the dual Higgs model with dual Dirac strings [11,12] the mass 
of a dual-vector field Mc is not proportional to the order parameter < XX > an d does 
not vanish in the limit < \X >— > 0. This is seen from the mass formula [4] 

M a (8M 2 +3 Ml) = -5Qg 2 < XX >, (4-2) 

which can be derived from Eq. (|1.2|) and the gap-equation Eq.( |1.4| ). Thus, in the MNJL 
model the dual-vector field does not need a Goldston boson as a longitudinal component. 
This distinguishes the transition to the non-perturbative superconducting phase in the 
NMJL and the dual Higgs model. Indeed, in the MNJL model this transition does not 
accompany the appearance of Goldston bosons. The former is rather natural, since the 
starting U(l) magnetic symmetry in the MNJL model is global and unbroken in the non- 
perturbative superconducting phase. Recall that in the dual Higgs model the magnetic 
U(l) symmetry is local and becomes spontaneously broken in the superconducting phase. 

Due to the independence of the mass of the dual-vector field on the monopole conden- 
sate the string tension o"o calculated in the tree-approximation does not depend on the 
monopole condensate too. The mass of the Higgs field M a replaced the cut-off Aj_, i.e. 
A± = M a . The dependence on the magnetic monopole condensate appears by virtue of 
the contributions of the quantum field fluctuations of the dual-vector C M and the scalar 
a fields around the shape of the Abrikosov flux line. Very similar to the dual Higgs 
model with dual Dirac strings the quantum field fluctuations increase the value of the 
string tension. This implies that for the consistent investigation of the superconducting 
mechanism of the quark confinement within the dynamics of magnetic monopoles and 
dual Dirac strings one cannot deal with a classical level only and quantum contributions 
should be taken into account. The string shape fluctuations of dual Dirac strings induce 
a Coulomb-like universal contribution calculated for openned strings by Liischer et al. [8] 
and for closed strings by Faber et al. [7]. 
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